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The influence of plasma on different effects of gravitational lensing is reviewed. Using the 
Hamiltonian approach for geometrical optics in a medium in the presence of gravity, an exact 
formula for the photon deflection angle by a black hole (or another body with a Schwarzschild 
metric) embedded in plasma with a spherically symmetric density distribution is derived. 
The deflection angle in this case is determined by the mutual combination of different factors: 
gravity, dispersion, and refraction. While the effects of deflection by the gravity in vacuum 
and the refractive deflection in a nonhomogeneous medium are well known, the new effect is 
that, in the case of a homogeneous plasma, in the absence of refractive deflection, the gravi¬ 
tational deflection differs from the vacuum deflection and depends on the photon frequency. 
In the presence of a plasma nonhomogeneity, the chromatic refractive deflection also occurs, 
so the presence of plasma always makes gravitational lensing chromatic. In particular, the 
presence of plasma leads to different angular positions of the same image if it is observed 
at different wavelengths. It is discussed in detail how to apply the presented formulas for 
the calculation of the deflection angle in different situations. Gravitational lensing in plasma 
beyond the weak deflection approximation is also considered. 


1. INTRODUCTION 

A gravitational lensing problem combines a well developed theory and a wide range of obser¬ 
vational phenomena connected with the deflection of light rays by gravity. The gravitational lens 
theory usually deals with a geometrical optics in vacuum and uses a notion of the deflection angle. 
Basic assumption is an approximation of a weak deflection angle of a photon. General Relativity 
predicts that a light ray passing near a spherical body of mass M with large impact parameter b is 
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Figure 1. Gravitational deflection of the light ray by the massive body, in vacuum. The light ray passing near the 
spherical body of mass M with impact parameter b is deflected by the small angle a, see (1). This picture 
corresponds situation when the impact parameter b is much larger than the Schwarzschild radius of the gravitating 
body. Einstein deflection does not depend on the photon frequency. 


deflected by the small angle (Fig. 1): 


a = 


2 Rg 
b 


4 M 

T’ 


G = c= 1. 


(1) 


This expression is valid if 6 > Rg, where Rg = 2 M is the Schwarzschild radius of the gravitating 
body. Deflection angle (1) is usually called ’Einstein angle’. In the most astrophysical situations 
related with gravitational lensing the approximation of weak deflection is well satisfied. Trajectories 
of photons in vacuum, as well as deflection angles, do not depend on the photon frequency (or 
energy), so gravitational lensing is vacuum is achromatic. 

The first gravitational lens outside the Solar system was discovered by Walsh et al [1], as two 
images of the quasar 0957 + 561. The correlated variability of these images had confirmed that 
both images belong to the same object, and are formed due to light deflection by the gravitation of a 
massive intermittent object. Tens of gravitational lens examples have been discovered in subsequent 
years, see, for example, [2]. 

In this paper we review an influence of plasma on different effects of gravitational lensing [3], [4], 
[5], [6]. The principal feature is that the deflection angles are changed in presence of plasma, and 
chromatic deflection occurs. In particular, it leads to angular difference in positions of images at 
different wave frequencies. For the light rays propagating through cosmic plasma such effects may 
be significant only for very long radiowaves. 1 

There exists another possible reason of curving of light trajectory: refractive deflection. In 


1 In vacuum the chromatic effects of gravitational lensing take place in study of microlensing of quasars by individual 
stars in lens (galaxy). Physical sources like quasars may have different sizes for different wavelengths. Magnification 
and the light curves expected as a result of microlensing depend on the source size, so chromatic effects arise. Such 
chromatic effects can give information about source structure, see [7], [8]. Sometimes the term ’chromatic gravi¬ 
tational lensing’ is used for this situation. Rigorously speaking, this situation is ’achromatic lensing of chromatic 
source’ rather than ’chromatic gravitational lensing’. In this paper we consider chromatic effects which related 
with the situation when the photon deflection itself becomes chromatic, in particular, the gravitational deflection 
of photon in plasma is chromatic. 
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non-homogeneous transparent media the photon moves along the curved trajectory, it is called as 
refraction. Refraction is widely known phenomenon, and it is easily observed at boundaries between 
water, air or glass, which is described by Snell’s law. Working in the frame of geometrical optics, we 
can characterize the properties of the medium by its refractive index. In inhomogeneous medium the 
refractive index depends explicitly on the space coordinates, and shape of trajectory is determined 
by this dependence. The refractive deflection of light has no relation to relativity and gravity, and 
takes place only in non-homogeneous media. If medium is dispersive too, which means that the 
refractive index depends on the frequency, then the refractive deflection is also chromatic. 

Considering propagation of a light under action of both gravity and plasma non-homogeneity, 
in linear approximation we can just summarize two deflection angles: a vacuum deflection due to 
gravitation of a point mass, and a refractive deflection due to the non-homogeneity of the plasma. 
These effects are separated from each other. This approach was used in consideration of deflection of 
radio rays passing the Sun [9] where the deflection was a separable combination of general relativistic 
deflection in vacuum and refraction in the coronal electron plasma, see also [10], [11]. The same 
approach is used in book of Bliokh and Minakov [12] where gravitational lensing was considered, 
produced by a gravitating body surrounded by a plasma with a spherically symmetric distribution. 

The approach above, considering propagation of light in presence of gravity and medium as 
a simple sum of two deflections, does not contain nonlinear effects. Possible mutual influence of 
medium and gravity is not taken into account. To consider this problem more rigourously, we need 
theory which includes gravity and medium presence simultaneously. 

A self-consistent approach to the light propagation in the gravitational field, in presence of a 
medium, was developed in a classical book of Synge [13]. For discussion and application of the 
Synge’s general relativistic Hamiltonian theory for the geometrical optics see also papers of Bicak 
and Hadrava [14], Kichenassamy and Krikorian [15], Krikorian [16]. 

Comprehensive review of general-relativistic ray optics in media is presented in the monograph 
of Perlick [17]. He used a different type of the approach to the problem. As one of applications of 
his formalism, general formulae for the light deflection angle in the Schwarzschild and Kerr metrics, 
in presence of spherically distributed plasma, have been obtained in the integral form. Formula for 
the Schwarzschild metric was later rederived on basis of Synge’s approach in the paper [6]. 

Investigations of non-linear effects connected with combined action of gravity and plasma were 
performed in papers of Bisnovatyi-Kogan and Tsupko [3], [4], [5], [6], on the basis of Synge’s theory. 
The main idea was to investigate the influence of plasma on the light trajectory in a homogeneous 
plasma where there is no refractive deflection. It was shown that the gravitational deflection itself 
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in the medium is not the same as in vacuum. Therefore, if we would like to consider tiny effects 
we should use corrected formula for the gravitational deflection instead of the vacuum gravitational 
deflection. 

This idea is the clearest in case when plasma is homogeneous. In this case the refractive deflection 
is absent. The question can arise: is the gravitational deflection of the light rays in this case the 
same as in vacuum? The answer is the following: the deflection will be differ from vacuum case [3]. 
For static plasma with the refractive index 



in a static weak gravitational field of a point mass, the gravitational deflection angle will be described 
by the simple formula [3]: 


s -t( 1 + i-|«)' (3) 

Here the frequency of the photon u>(r) depends on the space coordinate r due to the presence of a 
gravitational field (gravitational redshift). We denote cj(oo) = cu, e is the charge of the electron, m 
is the electron mass, oj e is the electron plasma frequency, Nq = const is the electron concentration 
of a homogeneous plasma. This formula is written for homogeneous plasma and valid under the 
condition of smallness of a (b 2> Rs )• At u e = 0 (concentration Nq = 0) or cj —>• oo this formula 
turns into the deflection angle (1) equal to 2Rg/b. 

Thus, it was shown for the first time in [3], [4] that due to dispersive properties of plasma even 
in the homogeneous plasma the gravitational deflection differs from vacuum deflection angle, and 
gravitational deflection angle in plasma depends on frequency of the photon. By using Synge’s 
approach, it was also found [3], [4] that the photon deflection is the same as in vacuum, if medium 
is homogeneous and non-dispersive (refraction index n=const). The physical reason is a dependence 
of the wave frequency on space coordinates in presence of gravity (gravitational redshift). In the 
paper [4] we have considered the gravitational lensing of radiowaves by the point or the spherical 
body in presence of a homogeneous and non-homogeneous plasma, together with the refractive 
deflection due to the plasma inhomogeneity, and comparison of these two effects have been there 
performed. All these results have been obtained in the approximation of smallness of a. 

Er and Mao [18] have studied the gravitational lensing when non-homogeneous plasma surrounds 
the lens. Using formulas for the deflection from [3], [4], they have calculated the extra deflection 
angle induced by the plasma for non-homogeneous plasma models. They have performed numerical 
modelling of different effects caused by this extra deflection. In particular they have shown that a 
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plasma presence may cause significant effects to be detected in low frequency radio observations (a 
few hundred MHz): the change of position due to a plasma can reach a few tens of milli-arcsec. 

Recently Morozova et al [19] have considered a gravitational lensing by a rotating massive object 
in a plasma. Their approach is based on papers [3], [4], they additionally take into account the 
rotation of lens and investigate new effects connected with presence of rotation. 

Theory of gravitational lensing usually considers small deflection angles (d <C 1), see Fig. 2. 
Light rays which are deflected by angles about 2n, 47T, ... can also form images. If the photon 
impact parameter is close to its critical value, a photon which goes from infinity can perform several 
turns around the central object and then go to infinity. Such photons form images (Fig. 3), which 
are called relativistic images [20]. Influence of plasma on properties of relativistic images was studied 
for the first time in [6]. 

The paper is organized as follows. In Section 2 we briefly describe the basic effects of gravitational 
lensing. In section 3 we consider the gravitational lensing in the case of a strong deflection of photons 
in vacuum. In section 4 we derive an exact formula for the deflection angle of photons in the 
Schwarzschild metric with spherically symmetric distribution of plasma. In section 5 we derive the 
exact formula for the angle of gravitational deflection of light in a homogeneous plasma. In section 6 
we consider the gravitational deflection in the homogeneous plasma in the case of a weak deflection. 
In section 7 we consider weak deflection of the photon in presence of gravity and non-homogeneous 
plasma. In section 8 we consider the gravitational deflection of light in the homogeneous plasma in 
strong deflection limit and discuss influence of plasma on relativistic images. In section 9 we make 
concluding remarks. 


2. BASIC EFFECTS OF GRAVITATIONAL LENSING 

In this section we make a short review of basic effects of gravitational lensing. For details, see, 
for example, books and reviews [21] - [28]. 

Gravitational lensing changes an apparent angular position of the source. Example of calculation 
of the trajectory of the photon for lensing with the small deflection is presented on Fig. 2. In this 
example the lens is a point body. Light ray from the source deflected by angle a by the point-mass 
gravitational lens with the Schwarzschild radius R$ goes to the observer. The observer sees the 
image of the source at angular position 0 which is different from the source position. Bending of light 
from the distant stars by gravity of the Sun leads to apparent displacement of the stars position. 
It can be observed for stars near Sun disc during eclipse, and it was one of the first observational 



6 



Figure 2. Change of apparent angular position of point source due to gravitational lensing. Trajectory of light ray 
is calculated in the Schwarzschild metric. Light ray from the source deflected by angle a by the point-mass 
gravitational lens with the Schwarzschild radius Rs goes to the observer. The observer sees the image of the source 
at angular position 9, which is different from the source position. R is the closest point of trajectory to gravitating 
center, it is usually referred as the distance of the closest approach, b is the impact parameter of the photon, Dd is 
the distance between observer and lens, D s is the distance between observer and source, Dds is the distance 
between lens and source. We should comment that at this picture the deflection angle is not small actually, and can 
not be calculated by approximate formula (1), we draw situation with such angle for visualization. 

proof of General Relativity [29]. 

Gravitational lensing leads to formation of multiple images, lensing on one gravitating point 
mass leads, in general, to formation of two images. On Fig. 4 there is the example of the simplest 
model of lens, the Schwarzschild point-mass lens. Due to bending of light the observer sees two 
images of the same source at different sides from lens. If source is point-like, then images of source 
are two points. If the source has a finite size, the deformation of images occurs, which is shown on 
Fig. 5 and Fig. 6 for source with circular shape. If source, lens and observer are aligned enough, the 
deformation can be rather big, and observer sees two arcs. Deformation becomes bigger, if angular 
position of source becomes closer to angular position of lens. In the case of perfect alignment, 
observer sees image in the form of ring - Einstein (Einstein-Chwolson) ring. If the gravitational 
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Figure 3. A scheme of formation of the relativistic image of the point source. Trajectory of light ray is calculated 
in the Schwarzschild metric. Before reaching observer, light ray perform turn around of black hole, at radius close 
to tm = SRs/'Z- Deflection angle of the photon in this case equals to d = 27r + Ad. There exist rays which make 
more revolutions around of black hole and form other relativistic images, the closest approach distances R for these 
rays are closer to tm- Combining with relativistic images at another side from the lens, these images form two 

infinite sequences of images. 

lens is not point-like, and has more complicated mass distribution, observer can see more than two 
images of the same source. In modelling of real system with multiple images the point mass lens is 
never used. Observed gravitational lens system are described usually by lens with broken spherical 
symmetry. 

Another effect of gravitational lensing is a so called magnification effect [21]. The surface bright¬ 
ness for the image is identical to that of the source in the absence of the lens. The flux of the image 
is the product of its surface brightness and the solid angle it subtends on the sky. The flux from 
source is changed due to gravitational lensing, because angular size of source is changed. Angular 
size (and flux) of each image is different from a source, and can be larger and smaller. The magni¬ 
fication factor fi, is the ratio of the flux of the image to the flux of the unlensed source, and is equal 







Figure 4. Lensing of the point source by the Schwarzschild point-mass lens in vacuum. Light rays from the source 
S are deflected by angle a by the point mass M. Observer O sees two images at opposite sides relatively to lens, 

instead of the real source S. The image which situated at the same side from the lens as real source is called 
primary image, another one is secondary. Rays from source which are not in plane of Source-Lens-Observer will not 
go to observer, therefore in case of point source the observer sees two images as point images. 

to the ratio of the solid angles of the image, and the unlensed source: /x = Au}/(Aoj)q. 

Gravitaional lensing leads also to time delay: in presence of lens the time of light travelling 
from source to observer is changed [21]. Time delay consists of two parts: geometrical delay, due 
to changing of trajectory shape, and potential or Shapiro delay, due to propagation of light in 
gravitational field. In the case when multiple images are observed, every image has its own delay, 
and time delay between images can be found. It is interesting as a way of finding of Hubble constant, 
see [30], [22], 

From the observational point of view there are three types of the gravitational lensing. 

(i) Strong lensing. In this case effect of lensing is strong what means the multiple imaging or 
strong deformation of image. Multiple images are close to each other at sky and have the same 
redshifts, what, together with spectral and time variation properties of the images, permits to obtain 
a clear conclusion about the presence of the gravitational lens. Extended arcs and rings are also 
examples of strong lensing. 

(ii) Weak lensing. In this case the lensing images are characterized by weak distortions and small 
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Figure 5. Primary and secondary images in case of lensing by the Schwarzschild point-mass lens in vacuum. 

Source is extended (not point but circle) and has internal structure drawn by curves inside circle, see picture. 

Observer does not see real source and sees two prolongated images with another angular sizes. At this picture the 

angular size of Einstein ring is also presented. 

magnifications. Usually only one image of each source is visible. In this case we can not identify 
fact of lensing using individual source because we do not know a real position of the source, its 
shape, and its real luminosity. To conclude the presence of lensing in this case is possible only in 
a statistical sense. Observations of many weakly distorted images of galaxies give a possibility to 
reconstruct the mass distribution in the gravitational lens. The example of the reconstruction of the 
distribution of the dark matter gravitational potential in merging Bullet cluster of galaxies 1E0657- 
558 (z~ 0.296) is given in [31]. It was obtained due to observations by HST/ASC of the field galaxies 
behind the clusters. This distribution is obtained together with the distribution of ordinary matter, 
by Chandra X-ray observations of the shocks formed in the merging cluster 1E0657-558, which 
brightness maxima do not coincide with the maxima of the dark matter gravitational potential. 
These observations are considered as a direct empirical proof of the existence of a dark matter. 

(iii) Microlensing. It is version of strong lensing in which the image separation is too small to be 
resolved. Lens which is a single star of about solar mass can lead to images split into microimages 
separated by microarcseconds (it explains an origin of the term). Multiple images due to stellar-mass 
objects can therefore not be resolved. The only observable effect of microlensing is the alteration of 
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Figure 6. Formation of arcs in case of lensing by the Schwarzschild point-mass lens in vacuum. Different angular 
positions of source relatively to lens are presented. Primary and secondary images become more prolongated if 
source is closer to lens (closer to line passing through observer and lens). In case of perfect alignment the observer 

sees ring. 

the apparent luminosity of background sources (change of flux), see [32], [33], [21]. There are several 
applications of microlensing effect. In the case of a quasar microlensing, multiple images (which 
occur due to strong lensing) of the quasar are studied. The images of a multiply imaged quasar are 
seen through a galaxy. Since galaxies contain stars, stellar-mass objects can affect the brightness of 
these images, and this is effect of microlensing for quasars. One of the most interesting application 
of microlensing is a microlensing by planets, which started with paper of Mao and Paczynski [34], 
for observation discovery see, for example [35]. 

From physical point of view we can distinguish situations of lensing on basis of magnitude of 
deflection angle: weak deflection of photon (small deflection angles, a <C 1) and strong deflection 
of photon (deflection angles are not small, a > 1). 

In most problems of the gravitational lensing the weak deflection approximation is applied. It 
means that the formula of Einstein deflection, in the case of a point lens, and formulae based on it, 
in the case of a mass distribution, are used. The lens equation [21] connects positions of source and 
images, at given geometry, in flat or expanded space-time. With a given mass distribution in the 
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Figure 7. Einstein ring, and relativistic rings. If the source, the lens, and the observer are perfectly alighned, a 
circle, known as the Einstein ring, is formed. Inside this ring there are rings formed by photons which have been 
deflected by angles about 2n, 47r, 67 t ...; these rings are sometimes referred to as relativistic rings. 

lens, the lens equation allows us to find positions, forms and magnifications of images, analytically 
or numerically. 


3. GRAVITATIONAL LENSING IN VACUUM IN THE CASE OF A STRONG DEFLECTION 

Weak deflection approximation is enough for the most problems and applications of gravitational 
lensing. The example of situation, when the deflection angles of photons beyond usual weak de¬ 
flection approximation are considered, is a lensing by black hole. In this case the photon can move 
close enough to central object to undergo large deflection. 

Let us consider the motion of a photon in the neighborhood of a black hole with a Schwarzschild 
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metric. The Schwarzschild metric is given by 

ds 2 = — (1 — \ dt 2 H-—---h r 2 (d6 2 + sin 2 6 dtp 2 ), Rs = 2 M, G = c = 1. (4) 

\ r / 1 — 2 M / r 

The shape of the orbit of a photon incident from infinity on a black hole is determined by its 
impact parameter b [36]: 

(i) if b < 3\/3 M, then the photon falls to Rs = 2 M and is absorbed by the black hole; 

(ii) if b > 3\/3 M, then the photon is deflected by an angle a and flies off to infinity. In 
gravitational lensing we are interested in this case. 

Exact expression for the deflection angle a in the case of a motion in the Schwarzschild metric 
can be derived from equations determining the photon orbit, see, for example, [36, 43]. For a given 
mass M, the deflection angle of a photon is a function of the radius of the closest approach R, and 
is represented by the integral 


a = 2 


dr 


R 


i JT _ Wi_ 2M) 
V b 2 r 2 \ r ) 


— 7r. 


The impact parameter 6, corresponding to the distance of the closest approach R is written as 


b 2 = 


R 3 


R - 2 M 


(5) 

( 6 ) 


The expression of the integral (5) in terms of elliptical integrals was first obtained in [37]. The 
expansions of exact deflection angle in powers of M/R and M/b is given in [38]. 

In the case of large impact parameters b 3> 3\/3Af the orbit is almost a straight line with a 
small deflection by an angle a = AM/b (see Fig.l). In this case (weak deflection limit) the impact 
parameter and the distance of closest approach are almost the same (see Fig.2). We can neglect a 
difference between the impact parameter and the distance of closest approach, and write formula (1) 
either with b or R. A form of the formula (1) in textbooks depends usually on a way of derivation 
(compare, for example, [39] and [36]). This case (weak deflection limit) is considered usually in the 
theory of gravitational lensing. 

There is another limiting case in which the photon deflection angle can be written in the analytical 
form. It is the case when the deflection angle is large, it is usually referred as strong deflection limit. 
If the value of the impact parameter is close to the critical value, 0 < b/M — 3\/3 <C 1, then the 
photon makes one or several turns around the black hole near a radius tm = 3M and flies off to 
infinity (Fig.3). In this case (the strong deflection limit) the deflection angle is [37], [40] 

a = —2 In ( —— 1 ] + 2 ln[12(2 - Vs)} - vr = 

W / 

R - 3 M 
36(2 - y/S)M 


7T = —2 In 


(i? — 3M)(2 + \/3) 


-2 In 


36M 


— 7T 


(7) 
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or, as a function of the impact parameter b [40], [41] 


a = — In y —— lj + ln[216(7 — 4-\/3)] — n = 

. b-3y/3M (b — 3\/3M)(7 + 4\/3) /o . 

648\/3(7 — 4\/3)M 648\/3Af 

In the strong deflection approximation, the deflection angle diverges logarithmically while the 
impact parameter approaches its critical value [40]. It was shown in [41] that for any spherically 
symmetric space-time the deflection angle diverges logarithmically when the impact parameter 
approaches its critical value. 

In the case of an unbound motion of massive particle near the black hole, the deflection angle is 
also written via integral and can be expressed via elliptic integrals (see, for example, [42]). In case 
of weak deflection (a< 1) the deflection angle is described by the formula (see, for example, [36]): 


Here b and v are the impact parameter and velocity of incident test particles, respectively. Analytical 
formulae for the case of a massive particle motion near Schwarzschild black hole, similar to (7), (8), 
have been obtained by Tsupko [42], in the limit of a strong deflection angle (a > 1), when the 
impact parameter approaches its critical value. 

Photons from a distant source which undergo one or several loops around the central object 
(lens), and then go to observer form images, which are called relativistic images [20] (see Fig.3). 
When the lens, source and the observer are situated on one line, these images form concentring rings 
around the lens [43], see Fig.7. In the opposite case the relativistic images are represented by two 
infinite sequences of ’spots’ situated on one line at both sides of the lens, and the closest approach 
distance of photons is converging to the radius 3 M. Using an exact expression for the deflection 
angle, Virbhadra and Ellis [20] calculated numerically the positions and the magnifications of the 
relativistic images for the Schwarzschild space-time. Frittelli, Kling, Newman [44] considered the 
exact lens equation for the Schwarzschild metric, and obtained solutions in the form of an integral 
expressions. Exact gravitational lens equation in spherically symmetric and static space-times was 
also considered in Perlick [45]. A review and comparison of different approximate lens equations is 
given in [46]. 

Relativistic images can be studied analytically with using of the deflection angle in the strong 
deflection limit. To make several revolutions around black hole, the distance of the closest approach 
must be very close to rM = 3 M, and formulae (7) and (8) provide very good accuracy. Bozza 
et al [40] applied analytical expression for the investigation of the relativistic images, finding the 
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positions and the magnifications of the relativistic images in the Schwarzschild metric. Relativistic 
rings (see Fig. 7) for a Schwarzschild black hole lens were considered in more details in [43]. 

Recent decade lensing beyond the weak deflection limit was studied in many papers, where the 
approaches described above were applied to different problems [20], [40], [45], [6], [43], [24], [41], [23], 
[46] - [58]. Relativistic images in vacuum for different types of metric [23], [24], or using alternative 
theories of gravitation [52] have been studied. 


4. GENERAL FORMULAE FOR THE PHOTON DEFLECTION ANGLE IN THE 
SCHWARZSCHILD METRIC IN PRESENCE OF PLASMA 


In this section we derive the exact expression for the deflection angle in the case of a photon 
motion in the Schwarzschild metric with a spherically symmetric distribution of plasma. We use 
spherical coordinates (r, 9, p) and the gravitational field is not supposed to be weak. Indices are 

i,k = 0,1,2,3; a, f3 = 1, 2,3 (r, 6, </?); signature {-,+,+,+}. (10) 


Let us consider a static space-time with the Schwarzschild metric (Rg = 2 M, G = c = 1): 

fj r 2 o M 

ds 2 = gikdx l dx k = — A(r) dt 2 + ^ ^ + r 2 ( dO 2 + sin 2 Odp 2 ) , A(r) = 1-. (11) 


Let us consider, in this gravitational field, a static inhomogeneous plasma with a refraction index 

u 2 o 4-7re 2 Af(?’) 


n = 1 — 


[cu(r)] 2 ’ 


w e = 


m 


( 12 ) 


Here N(r) is the electron concentration in plasma. 

The general relativistic geometrical optics in a curved space-time, in a dispersive medium with 
an angular isotropy of the refraction index, was developed by Synge [13], see also [14]. Let us assume 
that the space-time gik is given, and the refraction index n of an isotropic transparent medium is 
given by a scalar function of x a , and the photon frequency ui(x a ). The Synge’s method is based on 
the so called medium-equation: 


n 


2 


1 + PiP 1 
(pkV k f ' 


(13) 


This formula connects the components of the photon momentum 4-vector (the photon energy and 
3-vector of the photon momentum) and refraction index n. In this equation the metric g t k and the 
refraction index n are assumed to be known, a phase speed w = 1/n and the frequency uj(x a ) are 
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measured in the local rest frame, V k is 4-velocity of medium. There is also a relation for the photon 
frequency [13] 

Pi V i = -hw(x a ), (14) 


where h is the Planck constant. 

For a static medium in a static gravitational field, we have the medium-equation in the form 

[13]: 


n 


2 


1 + 


PiP 1 

(poV-9 00 ) 


(15) 


and 


PoV~9 00 = —p°V~9oo = -Tuv(x a ), 


(16) 


At infinity, in a flat space-time, we have: 

Po = -p° = -tu*j, (17) 

where ui = u( oo). 

In order to apply Hamiltonian method, we rewrite the medium-equation (15) in the form 

H(x\pi) = 0, (18) 


where 


H(x\pi) 



(19) 


Here we define the scalar function H(x l ,pi) depending on x l and pi. 

The trajectories of photons, in presence of a gravitational held, may be obtained from the 
variational principle [13] 


5 



= 0, 


( 20 ) 


with the restriction (18). The variational principle (20), with the restriction H{x l ,Pi) = 0, leads to 
the following system of Hamiltonian differential equations [13]: 


dx l 

~d\ 


8H 
dp i ’ 


dpi 

dX 


with the parameter A changing along the light trajectory. 


dH 

dx l 


( 21 ) 
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In the case of a plasma with the refractive index (12), the condition (18) can be reduced, with 
using of (12) and (16), to the form 


H(x\pi) = 


1 

2 L' 


g lk PiPk + u 2 e h 2 


= 0. 


From (21) we obtain the system of equations for the space components x a , p a : 


dx Q 

~dX 


= 9 aP Pp , 


( 22 ) 


(23) 


^ = ~ \ 9%PiPk -\h 2 (< 4) o . (24) 

Equation for the time component dpo/dX = 0 in the static metric (11) leads to po = const along the 
trajectory. The equation (16) at infinity has the form (17), so we find that the constant po equals 
to po = —Hu. 

Let us find the equation of the trajectory, and the photon deflection angle, for a motion in the 
equatorial plane 6 = 7r/2 of the metric (11). Then sin 2 0 = 1, and components of the metric have a 
form 

9rr = V7T, 9™ = A ( r ), 500 = ~ A (r), 5°° = (25) 

A(r) A (r) 

g w = r 2 , 9 W = y^ 9 ee = r 2 , 9™ = ^. (26) 

From the equation for 6 we have from (23) 


d6 ee dd Q 

Tx = 9 or Tx =p 


(27) 


It follows that for a motion in the plane 6 = 7t/ 2 we have po = p e = 0. From the equation (24) for 
p v it follows py, = const. Without a loss of generality we can assume that > 0. The equations 
for r and (p from (23) may be written as 


dr rr dp 

Tx = 9 Pr ' Tx =9 p - 


Substituting the components of metric we obtain: 

dp _ p v 


1 


(28) 


(29) 


dr r 2 p r A(r) 

The second multiplier in the right hand side of this equation can be expressed from the medium 
equation for plasma (22), which may be written in the form 


g rr p 2 r + g w p 2 v + g 00 p 2 0 + h 2 u 2 e {r) = A(r)p 2 + + h 2 u 2 e (r) = 0 . 


Po 


(30) 
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We have then 

p r A(r) = ±Jpl - A(r) + h 2 uj%(r)j . (31) 

Substituting (31) into (29), we obtain the equation for trajectory of photon: 

(32) 

Assume that a photon moves in such a way that its (^-coordinate increases. Then ’plus’ sign in (32) 
corresponds to the motion with the coordinate r also increasing, and the ’minus’ sign to the motion 
with r decreasing. 

For a photon which moves from the infinity to the distance of the closest approach R (minimal 
value of the coordinate r), and then to infinity, we have a change of angular coordinate as 
R oo 

gf d \ +/gf dr 2 (33) 

M~A(r) (^ + ft 2 cu2(r)) r T Jp 2 0 - A(r) + h 2 u%(r)^J 




= 2 


R r Jpo ~ A (r) (yf + h 2 u 2 e (r)} 


The motion along a straight line corresponds to the change of the angular coordinate A ip = n. 
Then the deflection angle may be written as 


a = 2 


R r \lPo ~ A{r) (^ + ^ 2 o;2(r)) 


The deflection angle (34) depends on a mass of the central body M, a distribution of plasma 
N(r), represented by w e (r), and on the parameters R, po and p ip . These parameters are connected 
by the boundary condition. The point r = R is a turning point, therefore in this point: dr/dX = 0 
and p r = 0. Using (31), we have in this point 

pI = A{R) (^ + hW e (R)^j . (35) 

Therefore only two parameters from {R, po, p^} are independent, while the third one is expressed 
through two others. The parameter p$ represents the photon energy at infinity. It is convenient to 
exclude p ip , and to derive the angle a as a function of po = —hw and R. We have from (35) 

2 _ d2 2 f 1 RiR) \ ('qtj’i 

P v ~ R Po[ TTm , .2 • 
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Substituting p 2 in (34) and using notation of Perlick [17] 


h(r) = r- 



uj 2 (r) 

1 / Mr) 

u 2 


I r cu 2 (r) 

i — 2 M uj 2 ’ 


we obtain the equation of the trajectory in the Schwarzschild space-time 

^ = ±- 1 (38) 

and the deflection angle of the photon moving from infinity to the central object and then to infinity 

OO 

a = 2 f - dr - -ir (39) 

i 

This expression for a was derived earlier in [17] by an another way. At given M and cu e (r), the 
deflection angle is determined by the closest approach distance R and the photon frequency at 
infinity oj. This formula allows us to calculate the deflection angle of the photon moving in the 
Schwarzschild metric in presence of a spherically symmetric distribution of plasma. 


a = 2 


5. GRAVITATIONAL DEFLECTION OF LIGHT IN HOMOGENEOUS PLASMA: EXACT 
EXPRESSION. ANALOGY BETWEEN A PHOTON IN A HOMOGENEOUS PLASMA 
AND A MASSIVE PARTICLE IN A VACUUM. 

Let us consider a homogeneous plasma, u e (r) = ui e = const. Rewriting the equation of the 
trajectory (32) and the expression for deflection angle (34), introducing notations of E and L 

^ = — = E > 1, ^ = L> 0, (40) 

nw e co e rvjje 

we obtain the equation of the trajectory as 

^ = ±4 , (41) 


/u 2 -A(r)( l + £) 


and the formula for deflection angle as 


a = 2 


R ^ sJE 2 -A{r)(l + £) 


The closest approach distance R and the parameters E and L are connected at the point r = R by 
the following boundary condition 

E 2 = A(R)( 1 + 4T) • (43) 
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The trajectory and the deflection angle are thus completely determined by any two parameters from 
{R, E, L}, with the third parameter being expressed through (43). Remind, that in vacuum the 
photon motion is determined by only one parameter, either R, or b, which are uniquely connected 
with each other, see (6). Transition to the vacuum from (43) can be done by tending E —>• oo and 
L —> oo (what corresponds to uj e = 0) and introducing the impact parameter as b 2 = L 2 /(E 2 — 1), 
as in the analogous case of the massive particle [36]. 

Plasma has unique dispersive properties, that leads to interesting analogy between motion of 
photon in homogeneous plasma and massive particle in vacuum. In the paper of Kulsrud and Loeb 
[59] (see also papers of Broderick and Blandford [60], [61]) it was shown that in the homogeneous 
plasma the photon wave packet moves like a particle with a velocity equal to the group velocity of 
the wave packet 



with a mass equal to the plasma frequency (multiplied by K) 

m v i = hw e , (45) 

and with an energy equal to the photon energy 

E p i = hw. (46) 

If we define E as an energy at infinity per unit rest mass, and L as the angular momentum per 

unit rest mass, then equations (41) and (42) describe the trajectory and the deflection angle of a 

massive particle with the rest mass huj e in the vacuum, see [36, 59, 62, 63]). 

As for massive particle in vacuum, bound elliptic orbits of the photon in homogeneous plasma are 
also possible. Discussion of gravitational binding of the photon in homogeneous plasma is discussed 
in [59],[6], 

The deflection angle (42) can be expressed via elliptic integrals. Due to described analogy, the 
expression for the deflection angle of a massive particle in vacuum is identical to (42), with an 
appropriate changes in the physical sense of E and L. The motion of a massive particles has been 
extensively studied, see [36, 37, 63-71]. In particular, the deflection angle of a massive particle 
can be expressed via elliptic integrals. Usually the deflection angle is rewritten as an integral of 
the polynomial of the third order in a radicand, see [6] for details. The form of the expression 
via elliptic integrals is determined by the roots of this polynomial. In a general case, these three 
roots are different, and supposed to be found numerically (or by cumbersome solution of cubic 
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equation, see, for example, [72]), for the given external parameters E and L. Depending on E 
and L, the roots have different relations between each other, and the expression of the deflection 
angle via elliptic integrals have therefore different form. In the newtonian limit these forms are 
expressed analytically, and are related to hyperbolic, parabolic, and elliptic motion of a massive 
particle around a gravitating center (Kepler problem), see [73]. The principal feature is that we 
cannot write the deflection angle as an explicit function of E and L. For given E and L, one has 
to obtain numerically the roots of the polynomial before being able to compute the deflection angle 
(see, for example, [67], [66]). 

In the paper [6] we suggested to use, as independent, another pair of parameters: R and E. 
Then one root of the polynomial is simply 1/R, and another two roots can be expressed analytically 
via R and E. This approach is convenient because we can perform all calculations analytically and 
obtain the deflection angle in the form of elliptic integrals as an explicit function of two parameters 
R and E. A similar approach is usually applied for the deflection of photons in vacuum, where an 
analogous polynomial is expressed via R. 

The problem is that R is not an ad hoc value and thus should not be chosen as an external 
parameter in applications. For example, for gravitational lens theory we need dependence of the 
deflection angle on the impact parameter. This difficulty can be avoided in strong deflection limit, 
where it becomes possible to express the deflection angle in terms of b and E, which is convenient 
for applications (see Section 8 and [6] for details). 

The deflection angle a for the photon in the homogeneous plasma can be written in the form [6] 


a = A x j — F(y,k) - ir. 


with y = 


8 MQ 


(6 M - R + Q)(R - 2M + Q) ’ 


(47) 


k = 


1 6M — R + Q 2 


2 Q 


, Q 1 = (R- 2 My + 8 M(R - 2 M) 


i , 2 M 

1 “T R(&- 1) 


Here F(z, k ) is an elliptic integral of the first kind [74, 75]: 

<p 


F{z, k) = j 


dx 


dQ 


^ ( l_ x 2 )(1 _ fe 2 x 2) J ^ 1 _ k 2 sin 2 0 
0 0 v 


= F(<P,k). 


(48) 


F(sin ip, k ) = F(ip, k) , x = sin 9 , z = sin ip. 
The expression a can be also written in another form [6] 


« = 4 t/ ^ i F 0-, F) - F(z, k)] - t r. 


(49) 
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Here F(l,k) = F(n/2,k) = K{k ) is the complete elliptic integral of the first kind, and 

2 _ 2 M + Q-R 
'“6 M+Q-R ' 


(50) 


This formula is written in the same form as the expression for the vacuum deflection angle, see 
[36], [64], [54], [43]. The difference between plasma and vacuum formulae is only in the expression 
for Q. When E goes to infinity, what corresponds to high energy photons, for which plasma effects 
are negligible, the expression for Q in (47) goes to Q 2 = (R — 2 M)(R + 6 M), and the formula (49) 
transforms into the formula for the vacuum deflection. The expression (47) for the vacuum case 
with Q 2 = (R — 2 M)(R + 6 M) is written as [37] 


a = 4\ 


■F(y,k) -7T, 


(51) 


with y = 


2 Q 


3 R -6 M + Q 


k = 


I6M-R + Q 


2 Q 


(52) 


The proposed novel technique can be also applied for a massive particle in vacuum, due to its 
analogy to a photon in plasma. In paper [42] we suggest to use either the pair (R, L ) or (R, E) as 
an independent pair of parameters, instead of (E, L). In that paper we obtain the deflection angle 
in the form of elliptic integrals as an explicit function of the two parameters (R, L) or (R, E). We 
present for the first time the analytical formulae for the deflection angles of massive particles in the 
strong deflection limit (a 2> 1) as an explicit function of external parameters at infinity, ( b , L) or 
( b , E). Parameters at infinity are simply related with each other, so it is possible to calculate the 
deflection angle with our formulae in very different situations. For example, our formulae allows one 
to calculate analytically the deflection angle of particle with given impact parameter and velocity 
at infinity as initial conditions, beyond usual small deflection case. 


6. GRAVITATIONAL LENSING IN THE LIMIT OF A SMALL DEFLECTION IN 
PRESENCE OF A HOMOGENEOUS PLASMA: ANALYTIC SOLUTION 

Exact expression for deflection angle in presence of Schwarzschild gravity and plasma (39) is 
manifestation of simultaneous and mutual action of different physical effects: gravity, refraction, 
dispersion. This angle (39) can not be calculated analytically. For analysis and applications of 
this result it is convenient to consider some particular cases with approximations. Considering 
case of small total deflection a <C 1, we have advantages to see clearly which physical reasons give 



22 


contributions to different parts of the total deflection. In particular, we will be able to discuss the 
difference of this angle and known, commonly used, formulae; see this and next sections. 

In the case of a homogeneous plasma the refractive index does not depend on space coordinates 
explicitly, so the refractive terms in deflection are absent. Therefore, in homogeneous case, we can 
consider this angle a as the gravitational deflection a gra v in given medium (plasma), and write 

® — OCgrav ( 53 ) 


We use notation a gra v here to emphasize that in homogeneous case the total deflection a consists 
of only the gravitational deflection (compare with the next section). 

Let us consider a situation when the closest approach distance is much larger than the 
Schwarzschild radius, R 2> M (Rs = 2 M). It means also that r 3> M (during the motion the 
r-coordinate changes from R to infinity), b 3> M, and the resulting deflection angle is small, a<l. 
Expanding the exact formula (39) with = const, we obtain the deflection angle in the form (see 
Appendix 1 for details): 


2 M (, 1 

a - a grav - — ^ + i _ w 2 /w 2 


or, in usual units, a = 


2GM 


1 + 


1 


Rs = 


2GM 


c 2 R \ \—(jj 2 /oj 2 / 

For gravitational lensing the dependence of angles on impact parameter b is needed, 
imation R M the difference between R and b is negligible, so we can just substitute 
(54), obtaining 


(54) 


In approx- 
R ~ b into 


a = 


2 M 
~b~ 


1 + 


1 — wf/w 2 


(55) 


The deflection angle (55) was directly derived in a simpler way in [3], [4], To do this, the 
Synge’s equations have been written in Cartesian coordinates. We have considered the photon with 
unperturbed trajectory in a homogeneous plasma as a straight line parallel to z-axis with an impact 
parameter b, and have found the deflection in the approximation of small perturbations, by the 
integration along the straight trajectory. Compare: the vacuum gravitational deflection in the form 
2 Rs/R was derived by expansion the integral of the exact deflection angle [39], and the deflection 
in the form 2 Rs/b was derived also by the integration of small deflections along the trajectory [36]. 

Formula (55) is valid only for u > oj e , because the waves with oj < u e do not propagate in 
the plasma. If oj ~ oj e , then gravitational deflection in plasma can be much larger than vacuum 
gravitational deflection, a 3> 2 Rs/b (condition a <C 1 must be satisfied too). 
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Formula (55) does not imply that u e /cj < 1. If it is supposed, we can rewrite the deflection (55) 


as: 


d — Qgrav — &einst &add — 


2 Rs + u 2 e 2 R s 


(56) 


b ' 2u 2 b ' 

Here we denote the vacuum gravitational deflection as a e i ns t and additional correction to the gravi¬ 
tational deflection connected with plasma presence as a ac id■ We emphasize that the correction a ac id 
is presented also in a non-homogeneous plasma, but has a more complicated form. 

The presence of plasma increases the gravitational deflection angle. In homogeneous plasma 
photons of smaller frequency, or larger wavelength, are deflected by a larger angle by the gravi¬ 
tating center. The effect of difference in the gravitational deflection angles is significant for longer 
wavelengths, when uj is approaching oj e . That is possible only for the radio waves. Therefore, the 
gravitational lens in plasma acts as a radiospectrometer [3]. 

We should use formula (3) when we consider gravitational lensing of radiowaves by point or 
spherical body in presence of homogeneous plasma. This effect has a general relativistic nature, in 
combination with the dispersive properties of plasma. We should also emphasize that the plasma 
is considered here like the medium with a given index of refraction, and this formula does not take 
into account gravitation of particles of plasma. 

The observational effect of the frequency dependence may be represented on the example of the 
Schwarzschild point-mass lens. Instead of two concentrated images with complicated spectra, we 
will have two expanded line images, formed by the photons with different frequencies, which are 
deflected by different angles (Fig.8). Estimations of an angular difference of positions of images in 
different bands can be found in our paper [4]. To see the difference in angular position of images, 
we should compare radio and optical observations of images. For optical frequencies the effect 
of plasma presence is negligible, and positions of images in this case can be calculated with the 
vacuum formulae. The difference between angular separations of images in vacuum and in plasma 
A0 q, produced by the same lens configuration with point-mass lens, is equal to 


A6» 0 1 w? 7 N e 

= -^~ 2 . 0 - 10 7 ^, 


Oo 


4 uj 2 


(57) 


where v is the photon frequency in Hz, ui = 27 tv, 9 o is angular radius of Einstein ring [4], This 
formula is written for homogeneous plasma, for non-homogeneous plasma see analytical formulae 
in next section of present paper and modelling and estimations in [18]. We can also compare two 
radio observations with the different frequencies. For one of two frequencies the effect of angular 
shifting is bigger than for another frequency, so we may observe the difference in positions of images 
of different frequencies. 
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Figure 8. Gravitational lensing in plasma, scheme of formation of the primary images of different frequencies. 

Two light rays which have different frequencies are shown. To reach the observer these rays are deflected by 
different angles. Therefore angular positions of images are different at different wavelengths. 

Gravitational lensing also leads to magnification of source. This means that the flux of image 
is bigger or smaller than the flux of source, different images have different magnifications. In 
observations we don’t know the intrinsic flux of source and its spectrum, but we can compare 
the fluxes of images in different frequencies. The ratios of the fluxes of images in different bands 
should be the same, if we consider lensing in vacuum, because gravitational deflection in vacuum is 
achromatic. 

The magnification is determined by the deflection law [21]. In the case of radio lensing in 
the plasma we have another formula for deflection instead of the Einstein angle, so formulae for 
magnification will be different. It leads to difference of magnigications of different images in different 
bands, when the light propagates in regions with different plasma density. By-turn it leads to 
difference of ratios of the fluxes of images in different bands (see details in [4]). It is another 
prediction of the model of gravitational radiospectrometer for observation. We should mention 
that Thompson scattering and absorption during propagation of radiation through the plasma can 
significantly change flux and complicate phenomena of lensing magnification. 
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Figure 9. Comparison of lensing by point source in vacuum (left) and in plasma (right). The photon trajectories 
and deflection angles in vacuum don’t depend on the photon frequency (or energy). Therefore lensing in vacuum is 
achromatic, and the Schwarzschild point-mass lens leads to two images with complicated spectra, formed by the 
photons of different frequencies which undergo deflection by the same angle. If lens is surrounded by plasma, the 
light rays of different frequencies have different deflection angles. The images become diffused and divide into many 
images formed by photons of different frequencies. Angular separation between images of different frequencies take 
place both in homogeneous and non-homogeneous plasma and is subject for observation. 

7. GRAVITATIONAL LENSING IN THE LIMIT OF SMALL DEFLECTION, IN PRESENCE 

OF A NON-HOMOGENEOUS PLASMA 

Discussion above is valid for the Schwarzschild lens in a homogeneous plasma. In the case of 
plasma non-homogeneity there is also the refractive deflection, a re f r . Our approach for gravitational 
lensing in plasma developed in [3], [4] allows us consider two effects simultaneously: the gravitational 
deflection in plasma which is different from the Einstein angle, and the refraction connected with 
the plasma inhomogeneity which does not depend on the gravity. In the paper [4] we have derived 
formulae for the deflection angle by the spherical body or spherically distributed gravitating matter 
in presence of a homogeneous and non-homogeneous plasma, with taking into account the gravitation 
deflection together with the refractive deflection. 

In applications the usual vacuum gravitational deflection a e i ns t gives the major contribution 
to the total deflection a, and all effects connected with plasma presence are significally smaller: 

add C^einst uud Ot re f r CTemst- 

Additional correction a a dd to the gravitational deflection caused by plasma (see previous section) 
is usually smaller than correction due to refraction a re f r , but there are situations when these cor¬ 
rections can be of the same order, see [4]. Let us consider situation when correction to gravitational 
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deflection is smaller than refractive deflection ( a a dd ay e /r)> and obtain the total deflection angle 
(in weak deflection approximation) as a sum of vacuum gravitational deflection a e i ns t and refrac¬ 
tive deflection a re f r . In this approximation we neglect a a dd and consider that the gravitational 
deflection is equal to vacuum gravitational deflection (a grav ~ oc e i ns t). 

Expanding formula for exact deflection (39), we obtain approximate formula in the case of a 
weak deflection (see Appendix 2 for details): 


Q! — Ot-einst &refr > 

(58) 

4 M 2Rs 

(59) 

&einst — ’ 


&refr 


oo 

RIU f 1 dN(r) , 

uj 2 J Vr 2 - R 2 dr 

R 


(60) 


For gravitational lensing the dependence of angles on the impact parameter b is needed. In 
papers [3] and [4] we have derive a using Cartesian coordinates. We have considered the photon 
with the unperturbed trajectory as a straight line parallel to z-axis, with the impact parameter b. 
We have obtained the gravitational deflection at given b as 

4 M 2Rg 

C^einst — ; — ; 

b b 



and refractive deflection as [3] 


or [4] 


K, 


d N 


a 


refr 


= -^ -77 dz 




db 


(62) 


OO 


K e b f 1 dN , 

&refr — n~ / 

cu z J r dr 


(63) 


o 

Formulae (60), (62), (63) are equivalent. In Appendix 2 we show how to transform formula (60) to 
(62) and (63). 

In the formula (62) N = N(r ) and r = Vb 2 + z 2 , so the expression under the integral sign is a 
function of b and z. To calculate the deflection angle, we perform partial differentation with respect 
to b at constant z and then perform integration with respect to z at constant b , and we obtain 
deflection angle as a function of b. 
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In the formula (63) we differentiate N = N(r) with respect to r, then substitute r = y/b 2 + z 2 . 
Expression under the integral sign will be a function of b and z , and we perform integration along 
z axis at constant b, and we obtain deflection angle as a function of b. 

All observational predictions about the difference in angular positions and fluxes in different 
frequencies take place also in the case of a non-homogeneous plasma distribution, because refractive 
deflection in plasma is also chromatic. We should emphasize that the presence of homogeneous or 
non-homogeneous plasma increases the gravitational deflection of photon. Vacuum gravitational 
deflection is usually considered as positive ( a e i ns t > 0 ) therefore the additional correction to grav¬ 
itational deflection due to plasma presence is also positive ( a a dd > 0). The refractive deflection 
can be both positive or negative, depending on the density profile. Usually the density of plasma 
in different models decreases with radius ( dN/dr < 0), therefore the refraction deflection is usu¬ 
ally opposite to the gravitational deflection: the correction due to refractive deflection is negative 
(a re f r < 0), see [4], 

For example, in the case of inhomogeneous plasma with a power-behaved concentration 

N(r ) = Vo(ro/r) fc , Nq = const, vq = const, k = const > 0 (64) 

the refractive deflection is ([3], see also [12], [9], [76], [11]) 



o 


Model with k = 1.25 is one of two models used in paper Er and Mao [18]. 

In the case of axisymmetric mass distribution the deflection with the impact parameter b is equal 
to the Einstein angle for the mass M(b), where M(b ) is the projected mass enclosed by the circle of 
the radius b [12, 21, 22, 77]. In another words it is mass inside the cylinder with the radius b. So, 
for this case, we should write: 


4 M(b) , 4 GM(b) 

oteinst = —7 - or, m usual units - 777 -. 

b c z b 


( 67 ) 
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8 . GRAVITATIONAL LENSING IN PLASMA IN THE LIMIT OF STRONG DEFLECTION: 

ANALYTIC FORMULAE 


Influence of plasma on the relativistic images have been considered for the first time in paper 
[6], where we have considered the case of a strong deflection angle for the light, traveling near 
the Schwarzschild black hole, surrounded by a uniform plasma. We derive for the first time the 
asymptotic analytic formulae for the gravitational deflection angle of photons in the homogeneous 
plasma, in a strong deflection limit. 

We have obtained the deflection angle a as a function of the closest approach distance R and 
the ratio of frequencies w 2 /w 2 in the form [6] 

a(R, x ) = —2 


1 + x 
2x 


In 


zi(x) 


R - r M 

rw 


— 7r. 


( 68 ) 


where 


zi(x) = 


9x — 1 + 2y / 6x(3x — 1) 


„, , 1 + x / 8 uSl 

, r M = 6M^-, x = ^l- — 


(69) 


48x 1 + 3x 

Here t'm is a critical (minimal) distance of the closest approach R, at given o; 2 /u; 2 . This formula 
(68) is asymptotic and is valid for R close to vm- At u 3 > id e the expression (68) becomes the 
vacuum formula (7). At the given ratio of frequencies cu 2 /u; 2 , the trajectory is determined by a 
choice of R. The deflection angle of a photon goes to infinity when R goes to rMi and photon 
performs infinite number of turns at circle with the radius t'm ■ If id 2> u> e , then x —)• 1, giving 
vm = 3 M, what corresponds to photons in the vacuum. The presence of a homogeneous plasma 
increases the radius of the critical orbits of the photons around a black hole. 

We have also obtained the deflection angle a as a function of the impact parameter b and the 
ratio of frequencies cu 2 /a; 2 : 


a(b, x) = — i 


,/L +I In 

2 zj(x) b — b cr 


V 2x 

3x b cr 



1 + x 


(70) 


3x- 1 

This formula is valid for b close to 6 cr , b cr is a critical value of impact parameter at given t u 2 /w 2 . 
For u 3> tu e we obtain the critical impact parameter for vacuum: b cr = 3\/3 M. At u > w e the 
expression (70) becomes the vacuum formula (8). 

We apply these formulae for the calculation of positions and magnifications of relativistic im¬ 
ages in the homogeneous plasma [6], which are formed by the photons performing one or several 
revolutions around the central object (see Fig. 10). For simplicity, let us rewrite formula (70) as: 

b b cr 


a(b,x ) = —a(x) In 


c(x ), 


(71) 
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where a(x) and c(x) are written as 


a ( x ) = \l-4—i c{x) = In 


2^(x) 

3x 


~ 7T, 


2x ’ ' 7 V 2.x 

zi(x), vm and x are defined in (69), 6 cr is defined in (70). 

To obtain the impact parameters corresponding to the relativistic images we can write an equa¬ 
tion: 


a(b,x) = 2nn, n = 1,2,... 


(72) 


Here n is number of pair of relativistic images. More rigorously, the deflection angle in this case is 
d = 27T +Ad (see Figures 3 and 10), where Ad has different values for different positions of source. 
But since Ad <C 27rn we can write equation for b in the form (72), with high accuracy. We obtain 
the impact parameters b n {x) of the relativistic images: 

c(x) — 27rn' 


bn — b r 


1 + exp 


a(x) 


The corresponding angular positions of the relativistic images are: 


a — cr 

0n -w d 


1 + exp 


c(x) — 2irn 
a(x) 


(73) 


(74) 


where Dj is the distance between observer and lens. 

The angular positions 9 n in plasma is always bigger than the positions in vacuum. In Table 1 
results for E = 2 (cj = 2co e ) are presented. Therefore the presence of homogeneous plasma increases 
the angular separation of the point relativistic images from gravitating center or the angular size of 
the relativistic rings. 

Let us consider b n near the vacuum values. With E 2 —>• oo (lo 3> u e ) we have: 


r M ^ 3 M ^1 + , b cr ~ 3\/3Af (l + , 0„ = , 

b n ~ 3\/3M(l + C c ) + 3^3{ 3 + C C [^(2n + 1) + 2^3 - 3] } , 


where 


C c = 216(7 - 4\/3) e ~ 7r ( 2n + 1 ) 


(75) 

(76) 

(77) 


Note that E = 00 corresponds to the vacuum. 

Let us consider a behavior of b n in opposite situation: when the photon frequency is close to 
plasma frequency. With E 2 —>• 1 (w ~ u e ) we have: 

AM AM 

_ r~^J _ 

VE 2 - 1 ripi 


r M — AM[1 + (E 2 — 1)] , b c 


here n p i = \ 1 - 






(78) 
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Figure 10. Gravitational lensing in homogeneous plasma. Scheme of formation of the relativistic images of 
different frequencies. The photons of smaller frequency are deflected by a larger angle by the gravitating center. 

The first ray has the frequency in this case plasma effects are negligible, and the trajectory can be 

computed using the vacuum equations (see Fig.3, geometry and positions of objects are the same). The other ray 
has the frequency oj ~ ui e , in this case plasma effects are significant, and the trajectory should be computed using 

plasma equations. 


9n = ^~, b n ~ ^ [l + 32 e -*(2n+i)/V2-] (79) 

L)d VE 2 — 1 L j 

We see that with uj —> oj e angular sizes 9 n increases unboundedly. 

Let us now consider the magnification of the relativistic images. We consider the relativistic 

images of a point source located at the angular position f3 from line connecting the observer and the 

gravitational center (lens). We use approximate approach, similar to [41], based on our analytical 

formula (71), see [6] for details. More precise way is to use the exact formula (39) and to solve this 

problem by the similar approach as it was done in [20] for relativistic images in vacuum. 

We obtain: 
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Table 1. Comparison of the impact parameters and the angular sizes of the relativistic rings in vacuum and 
homogeneous plasma for oj = 2ui e (E = 2). For convenience, differences b n — b C r are presented. The values of impact 
parameters are given in units of mass M. The angular sizes can be calculated as 9 = b/Dd . These values can also 
be considered as positions of two relativistic images of the point source, in case if there is no perfect alignment. 


Relativistic rings 

vacuum 

plasma 

Critical value ( b cr ) 

3 ^ 

3>/3 + 0.540300928156 

Ring 1 (bi - b cr ) 

0.00650 

0.00899 

Ring 2 ( b 2 - b cr ) 

1.21- 10~ 5 

2.05- 10~ 5 

Ring 3 (6 3 - b cr ) 

2.27- 10 -8 

4.68- 10 -8 


and a(x) and c(x ) are coefficients defined in (71). In expression (80) the variables b cr , l n , a(x), 
c(x) depend on x, so these variables depend on the ratio of the photon and the plasma frequencies. 
Magnification /u n of the relativistic images of the point source formally tends to infinity if the source 
angular position f3 goes to 0, as it takes place in vacuum [20]. 

In case of vacuum (ui/u} e = oo): 




vac 

n 


D s b 2 cr (l + r n ac )l h n 

D ds D 2 (3 


b cr = 3V3M , 


where 1™ C is defined in (77). 

If E 2 —> 1 (uj ~ u e ). we have: 


(81) 


D s 16 M 2 (1 + l P X 
>ln D ds D 2 E 2 -l ^2/3 


Ipi _ 32 e ~ 7r ( 2r! + 1 )/^ 


(82) 


We see that if u> —> u ei magnifications fi n can increase unboundedly. 

Numerical examples are given in Table 2. 

Possibility and difficulties of observation of relativistic images in vacuum are discussed in big 
details in [20] and [49]. Presence of plasma formally may lead to much stronger magnification of 
relativistic images at oj —>• u e , but increasing absorbtion creates additional problems for observations. 

While there are many observations of multiple ’ordinary’ images formed by photons weak 
deflected by massive lenses, the observations of relativistic images formed by the photons in strong 
deflection regime seems to be very unlikely. This low probability becomes even less, if we take into 
account that there are reliable observations of black holes only with accretion disc environment 
which significantly absorbes the radiation. 

Nevertheless, existence of isolated black holes in edge of Galaxy can not be fully excluded. 
Detection of them by usual methods is almost impossible, but lensing of the distant source radiation 
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(for example, star of Magellan Clouds) by such black hole might be observed with using some non¬ 
trivial technique. This detection is unlikely, but if it occurs, it would be very informative both for 
finding the position of lensing object and for properties of the surrounding medium. 

Table 2. Comparison of the magnification factors of relativistic images for lensing in homogeneous plasma and in 
vacuum. Ratios y,n/ y00 c are presented, at given D s , Dd, Dds, 0- Values and y00 c are calculated with (80) and 
(81) correspondingly. Values g n and g00 c depend on D s , Dd, Dd s , 0, u;/w e , but ratio g n /iJ00 c depends only on 
w/w e . In last column the values g00 c are calculated for M/Dd = 2.26467 • 10 11 , what corresponds to supermassive 
black hole in center of Milky Way, D a /Dds = 2, 0 = lgas (these parameters have been taken from [49]). 


Dn/t C° 

U>/oj e = 1.1 

w/w e = 2 

u>/u e — 10 

..vac 

rn 

Ml lkT° 

13.3 

1.48 

1.01 

0.716- 10" 11 

M2/Mr c 

40.0 

1.81 

1.02 

0.134- 10 -13 

M3 IpT 

120 

2.21 

1.03 

0.249- 10" 16 


9. CONCLUSIONS 

(i) The model of gravitational lensing in plasma is developed in details, by using Synge’s hamil- 
tonian theory for geometrical optic in gravity and medium. Exact formula for the deflection angle 
in the Schwarzschild metric, surrounded by plasma with spherically symmetric density distribution 
is written as (39), without the limit for the value of the deflection angle, the distance of the closest 
approach, or the ratio of plasma and photon frequencies. In the case of a homogeneous plasma the 
integral for the deflection angle is expressed via elliptic integrals, see Section 5. 

(ii) In presence of both gravity and plasma the deflection angle is physically defined by mutual 
combination of different phenomena: gravity, dispersion, refraction. Effects of deflection by gravity 
in vacuum and the refractive deflection in non-homogeneous medium are well known and presented 
in literature with necessary details. New effect, which we emphasize, is that in the case of a 
homogeneous plasma, in absence of refractive deflection, the gravitational deflection differs from 
the vacuum deflection, see (54) and (55). 

(iii) Plasma is a dispersive medium, its refractive index depends on the photon frequency. There¬ 
fore the deflection angle in presence of both gravity and plasma always depends on the photon fre¬ 
quency. In particular, presence of homogeneous plasma leads to chromatic gravitational deflection; 
the presence of non-homogeneous plasma leads also to chromatic gravitational deflection and to 
chromatic refractive deflection. So presence of plasm.a always makes gravitational lensing chromatic 
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and leads to difference in angular position of the sam.e image at different wavelengths. 

(iv) For applications to gravitational lensing problems connected with observations, see section 6 
and 7 where case of weak deflection is discussed. For the deflection by point-mass in a homogeneous 
plasma see (55), for the inhomogeneous case see (58) with (61), (60), (62), (63). In the case of 
gravitating mass distribution (in our formulae we neglect the mass of plasma particles) the projected 
mass M(b ) should be used (67), as it is usually done in lensing. As an example of application of 
these formulae, see numerical modelling of gravitational lens with plasma which is performed in [18], 
with attention to many gravitational lens phenomena and discussion of observational possibilities. 

(v) The deflection angles for light rays near Schwarzschild black hole in homogeneous plasma, 
in strong deflection limit, are derived. Using of strong deflection limit allows us to investigate 
analytically properties of relativistic images in presence of plasma. 
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APPENDIX A: APPENDIX 1 


Let us rewrite formula (39) as 


with 



dr 





— 7T , 


h{r) = r. 


/ i 

wf(r) 

i/ 1 - 2 M/r 

cu 2 


(Al) 


(A2) 


We consider case with = const. We denote: w 2 /cu 2 = B. Using that Ri$> M and r 3> M, we 
obtain: 


h 2 (r) ~ r 2 (l + ^-B) _ r2 ( 1 + r(!A)) 

mW) ~ W(\ + 'if - b) ~ R 2 (j + 




(A3) 
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1 + 


Transforming 


and 


2 M 


R 2 V r(l — B) R{l-B)J R 2 i? 2 1 - B ri? 


T-iO T 


r 2 2M R — r 


h 2 (r) 

-1/2 

h 2 (R) 

— 


R 


r 2 — R 2 | 2M r(i? — r) 
i? 2 + 1 - B R? 

M r(r — R) 


-1/2 


1 + 


Vr 2 - R 2 L l-B R(r 2 -R 2 )\ 


(A4) 


we get 


1 ( M 

— - ( 1 H- 


’\J 1 — 2M/r r \ r 


a = 2 


R 


dr + 2 


R 


ry/r 2 — R 2 J rVr 2 — i? 2 f 


oo oo 

f R M , f M r-R 

dr+ 2 I --- —_ _ , n dr -it 


R 


R 


1 _ B (r 2 - i? 2 ) 3 / 2 


(AS) 


(A6) 


2 M M 1 2M 1 

_ ^ + “r - + 2 i-s i~ 71 ~ v 1 + i-w 2 /w 2 / ■ 


APPENDIX B: APPENDIX 2 


Let us expand formula (39) with M/r <C 1 and w 2 (r)/w 2 <C 1. We denote w 2 (r)/w 2 = £> r and 
cj 2 (i?)/ur = Br. We get: 


OO OO 


a = ^straight + a einst + a re / r - 7T = 2 J fodr + 2 J j\ dr + 2 j f 2 dr -it , 


(Bl) 


where 


fo = 


R R 


R 


R 


fl = 


ry/r 2 — R 2 ’ 

(r 2 + rR+ R 2 )M 
r 2 (r + R)\fr 2 ^-R 2 ' 


(B2) 

(B3) 


h = 2(r 2 -^2)3/2 - B «>- < B4 > 

Integration of /o leads to a s t ra ight — 11 which is cancelled by ir in (Bl). Integration of f± leads to 
o^einst = 2 M/R which can be written here as 2 M/b. Integration of f 2 is: 

OO OO 

a-refr = J ^ 2 _ r2 ^ 3/2 ( B r ~ b r) dr = -R J (B r - B R ) d (^===j = (B5) 

R R 
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R(B r — Br) 


Vr 2 - R 2 


+ R 


R 


R 


1 dB r 
Vr 2 - R 2 dr 


dr. 


The first term is evidently equals to zero at r = oo. At r —>• R we can write: 


R(B r — Br) 


Vr 2 - R 2 


and we obtain: 


>R V 2R ( r ~ R ) 


R ( Br + 


dB r 


dr 


(r-R)- Br 


r=R 


oc \Jr — R 


r^rR 


r^fR 


= 0 , 

(B6) 


^refr R 


1 dB r RK e 
dr = —rr- 


\Jr 2 - R 2 dr 


ur 


1 dN(r) 
Vr 2 - R 2 dr 


dr. 


(B7) 


R R 

Let us transform this formula to variables b and z. We consider an unperturbed trajectory as 
straight line parallel axis z with impact parameter b , and r = \/b 2 + z 2 . We transform r-derivative 
to derivative with respect to perpendicular to z direction which we can denote as b: 


dN dN dr dN b 
8b dr db dr r 

We transform integration with respect to r to integration with respect to z, at constant b: 

z dz 


(B8) 


dr = d(\/b 2 + z 2 ) = 


(B9) 


We can change R to b. 
Substituting 


we obtain 


D I r~2 dN r dN zdz 

R = b, Vr 2 - R 2 = z , — = - —, dr = -, 

dr b db r 


bK P r lrdN z 


&refr ^2 


b db 


dz = 


oo 

K e r dN 




db 


dz 


(BIO) 


(BH) 


and 


OO 

K P b r 1 dN 


refr 9 / j 

u r r dr 


dz. 


(B12) 


1. D. Walsh, R. F. Carswell and R. J. Weymann, Nature, 279, 381 (1979). 

2. C.S. Kochanek, E.E. Falco, C. Impey, J. Lehar, B. McLeod, H.-W. Rix, CASTLES Survey, 
http: //www.cfa.harvard.edu/castles/ 



36 


3. G. S. Bisnovatyi-Kogan and O. Yu. Tsupko, Gravitation and Cosmology, 15, 20 (2009). 

4. G. S. Bisnovatyi-Kogan and O. Yu. Tsupko, Mon. Not. R. Astron. Soc. 404, 1790 (2010). 

5. O. Yu. Tsupko and G. S. Bisnovatyi-Kogan, Gravitation and Cosmology, 18, 117 (2012). 

6. O. Yu. Tsupko and G. S. Bisnovatyi-Kogan, Physical Review D 87, 124009 (2013). 

7. R. Kayser, S. Refsdal, and R. Stabell, Astronomy and Astrophysics, 166, 36 (1986) 

8. J. Wambsganss and B. Paczynski, Astronomical Journal, 102, 864 (1991) 

9. D. O. Muhleman, R. D. Ekers, and E. B. Fomalont, Phys. Rev. Lett. 24, 24, 1377 (1970). 

10. D. O. Muhleman and I. D. Johnston, Phys. Rev. Lett. 17, 8, 455 (1966). 

11. A. P. Lightman, W. H. Press, R. H. Price, and S. A. Teukolsky, Problem Book in Relativity and 
Gravitation. Princeton University Press, Princeton, New Jersey, New York, 1979. 

12. P. V. Bliokh and A. A. Minakov, Gravitational Lenses. Naukova Durnka, Kiev, 1989. (in Russian) 

13. J. L. Synge, Relativity: the General Theory. North-Holland Publishing Company, Amsterdam, 1960. 

14. J. Bicak and P. Hadrava, Astron. and Astrophys. 44, 389 (1975). 

15. S. Kichenassamy and R. A. Krikorian, Physical Review D 32, 1866 (1985). 

16. R. A. Krikorian, Astrophysics 42 (3), 338 (1999) 

17. V. Perlick, Ray Optics, Fermats Principle, and Applications to General Relativity. Springer-Verlag, 
Berlin, 2000. 

18. Xinzhong Er and Shude Mao, Monthly Notices of the Royal Astronomical Society, 437, 2180 (2014). 

19. V. S. Morozova, B. J. Ahmedov, and A. A. Tursunov, Astrophysics and Space Science, 346, 513 (2013). 

20. K. S. Virbhadra and G. F. R. Ellis, Phys. Rev. D 62, 084003 (2000). 

21. P. Schneider, J. Ehlers, and E. E. Falco, Gravitational Lenses. Springer-Verlag, Berlin, 1992. 

22. P. Schneider, C. S. Kochanek, and J. Wambsganss, Gravitational Lensing: Strong, Weak and Micro, 
Swiss Society for Astrophysics and Astronomy Series: Saas-Fee Advanced Courses, Number 33. Springer, 
Berlin, 2006. 

23. V. Perlick, Living Rev. Relativ. 7, 9 (2004) 

24. V. Perlick, eprint arXiv: 1010.3416 (2010) 

25. M. Bartelmann, P. Schneider, Physics Reports, 340, 291 (2001) 

26. M. Bartelmann, Classical and Quantum Gravity, 27, Issue 23, id. 233001 (2010). 

27. J. Wambsganss, Living Reviews in Relativity, 1, 12 (1998) 

28. H. Hoekstra, Bh. Jain, Annual Review of Nuclear and Particle Systems, 58, Issue 1, 99 (2008) 

29. F. W. Dyson, A. S. Eddington and C. Davidson, Philosophical Transactions of the Royal Society of 
London. Series A, 220, 291 (1920). 

30. S. Refsdal, Monthly Notices of the Royal Astronomical Society, 128, 307 (1964) 

31. D. Clowe et ai, The Astrophysical Journal, 648, Issue 2, L109-L113 (2006) 

32. S. Refsdal, Monthly Notices of the Royal Astronomical Society, 128, 295 (1964) 

33. A. V. Byalko, Astron. Zh. 46, 998 (1969), (English translation: Sov. Astron. 13, 784, (1970)) 

34. S. Mao, B. Paczynski, ApJ, 374, L37 (1991). 



37 


35. J.-P. Beaulieu, D. P. Bennett, P. Fouque et al , Nature 439, 437 (2006). 

36. C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation. Freeman, New York, 1973. 

37. C. Darwin, Proceedings of the Royal Society of London, Series A, Mathematical and Physical Sciences, 
249 (1257), 180 (1959) 

38. C. R. Keeton and A. O. Petters, Phys. Rev. D 72, 104006 (2005) 

39. L. D. Landau and E. M. Lifsliitz , The Classical Theory of Fields. Pergamon, Oxford, 1993. 

40. V. Bozza, S. Capozziello, G. Iovane, and G. Scarpetta, Gen. Relativ. Gravit. 33, 1535 (2001). 

41. V. Bozza, Phys. Rev. D 66, 103001 (2002). 

42. O. Yu. Tsupko, Physical Review D 89, 084075 (2014) 

43. G. S. Bisnovatyi-Kogan and O. Yu. Tsupko, Astrophysics, 51, 99 (2008). 

44. S. Frittelli, T. P. Kling, and E. T. Newman, Phys. Rev. D 61, 064021 (2000). 

45. V. Perlick, Phys. Rev. D 69, 064017 (2004) 

46. V. Bozza, Phys. Rev. D 78, 103005 (2008). 

47. V. Perlick, in Proceedings of the MG11 Meeting on General Relativity, edited by Hagen Kleinert, Robert 
T. Jantzen, editor of the Marcel Grossmann Meeting series: Remo Ruffini, published by World Scientific 
Publishing Co. Pte. Ltd., 2008, pp. 680-699, eprint arXiv:0708.0178 (2007) 

48. W. Hasse and V. Perlick, Gen Relativ Gravit 34, 415 (2002) 

49. K. S. Virbhadra, Phys. Rev. D 79, 083004 (2009). 

50. Clarissa-Marie Claudel, K. S. Virbhadra, and G. F. R. Ellis, Journal of Mathematical Physics 42, 818 

( 2001 ). 

51. K. S. Virbhadra, C. R. Keeton, Phys. Rev. D 77, 124014 (2008). 

52. V. Bozza, Gen. Relativ. Gravit. 42, 2269 (2010). 

53. E. F. Eiroa, G. E. Romero, and D. F.Torres, Physical Review D 66, 024010 (2002) 

54. S. V. Iyer and A. O. Petters, Gen Relativ Gravit 39, 1563 (2007) 

55. E. F. Eiroa and C. M Sendra, Class. Quantum Grav. 28, 085008 (2011) 

56. N. Mukherjee and A. S. Majumdar, Gravitation and Cosmology, 15 (3), 263 (2009) 

57. T. Ghosh and S. SenGupta, Physical Review D, 81 (4), 044013 (2010) 

58. Shao-Wen Wei, Yu-Xiao Liu, Chun-E Fu, and Ke Yang, eprint arXiv: 1104.0776 (2011) 

59. R. Kulsrud and A. Loeb, Physical Review D, 45, 525 (1992) 

60. A. Broderick and R. Blandford, MNRAS, 342, 1280 (2003) 

61. A. Broderick and R. Blandford, Ap&SS, 288, 161 (2003) 

62. S. Weinberg, Gravitation and cosmology: principles and applications of the general theory of relativity. 
John Wiley and Sons, Inc., New York London - Sydney - Toronto, 1972. 

63. Ya. B. Zel’dovich and I. D. Novikov, The Theory of Gravitation and Stellar Evolution. Nauka, Moscow, 
1971. (in Russian) 

64. S. Chandrasekhar, The mathematical theory of black holes. Clarendon Press Oxford, Oxford University 
Press, New York, 1983. 



38 


65. Yu. Hagihara, Japanese Journal of Astronomy and Geophysics, 8, 67 (1931) 

66 . A. W. K. Metzner, Journal of Mathematical Physics, 4 (9), 1194 (1963) 

67. B. Mielnik and J. Plebanski, Acta Phys. Polon. 21, 239 (1962) 

68 . C. Darwin, Proceedings of the Royal Society of London, Series A, Mathematical and Physical Sciences, 
263 (1312), 39 (1961) 

69. A. F. Bogorodsky, Einstein’s Field Equations and Their Application to Astronomy. Kiev Univ., Kiev, 
1962. (in Russian) 

70. C. Miro Rodriguez, Nuovo Cimento B, 98 (1), 87 (1987) 

71. C. Miro Rodriguez, Nuovo Cimento B, 100 (6), 801 (1987) 

72. I. M. Bronstein and K. A. Semendyaev, Handbook of Mathematics. Van Nostrand, New York, 1985. 

73. L. D. Landau and E. M. Lifsliitz, Mechanics. Nauka, Moscow, 1988. (in Russian) 

74. G. A. Korn and T. M. Korn, Mathematical handbook for scientists and engineers: definitions, theorems, 
and formulas for reference and review. Courier Dover Publications, 2000. 

75. I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, series and products. Academic Press, New York, 
1965. 

76. C. Thompson, R. D. Blandford, Cli. R. Evans, and E. S. Phinney, ApJ, 422, 304, (1994), Appendix A 

77. E. E. Clark, MNRAS, 158, 233 (1972) 



